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Abstract 

We study the melting of a moving vortex lattice through numerical simula- 
tions with the current driven 3D XY model with disorder. We find that there 
is a first-order phase transition even for large disorder when the correspond- 
ing equilibrium transition is continuous. The low temperature phase is an 

anisotropic moving glass. 

PACS numbers: 74.60.Ge, 74.60.Ec, 74.50+r 



1 



Typeset using REVTpjX 



There is now growing experimental [0-0] and theoretical [|5|-|TTH evidence that the vortex 
lattice (VL) in clean high T c superconductors melts via a first order transition. Most of 
the experiments consist of measurements of jumps in the linear resistance ]1]-[J and in the 
equilibrium magnetization Q . These results are consistent with computer simulations in the 
3D XY model p-pH[] as well as in other models flLl|| . However, even in clean samples there 



is always random weak pinning, which destroys the crystalline order of the VL. Moreover, 
large disorder can transform the first-order transition into a continuous transition, as found 
both experimentally || and in simulations ||. Koshelev and Vinokur []12 have proposed 



that melting from a perfect lattice will only be possible for a rapidly moving VL. For driving 
currents I much larger than the critical current I c , the VL is depinned and the effect of the 



in 



random potential in the moving VL is considerably weakened. Recent experiments ||13Hi5 
clean YBa2Cu307 crystals driven at currents I 3> I c have indeed shown a melting transition 
of the moving VL. 



On the other hand, Giamarchi and Le Doussal [16| have shown that the perturbation 
theory of [|T^] breaks down even at large vortex velocities, because some modes of static 
disorder are still present in the moving system. Balents and Fisher [ I7| , in the related 
system of sliding charge density waves, found that the moving phase has only quasi-long 
range order in 3D. Therefore, questions such us the existence and nature of a solid moving 
phase and which effects of the static disorder remain once the VL is in motion are currently 
under active discussion. 

In this Letter we study the behavior of a moving VL driven by currents / > I c as a 
function of temperature and disorder. We perform numerical simulations in a 3D current- 
driven XY model MM. A distinct advantage of this model with respect to 2D molecular 



dynamics simulations |r8|JT9[] , is that, besides the higher dimensionality, the model allows 
for local fluctuations of the vorticity (i.e. thermally nucleated vortex loops). Here, we find 
that the moving VL is a solid phase with an anisotropic structure factor, consistent with 



the theory of which melts via a first-order phase transition. 



The 3D XY Hamiltonian, 



^=-E J ACOs[A^(r)-A A (r)] , (1) 

r,jX 

results from modeling the macroscopic properties of a superconductor with the thermal 
fluctuations of the phase 6 [|BHT0|]. We consider a network || given by r = r\S\ + r 2 §2 + r 3 z, 
with §i = x, §2 = — ^x+^y, S3 = \x+^-y. The phase difference is A^6>(r) = 9(r+jl) — 9(r). 
The gauge factor A A (r) = ^ J r r+M A ■ <il ($ = h/2e) depends on the magnetic induction 
B = V x A = Bz. (We neglect fluctuations in B, valid for ($ /-B) 1//2 < A). We have 
A z (r) = and X^A^r) = 2nf = 27t_Bs 2 a/3/4$o, with the network discretization s = 
We model disorder with a random J A in the interval [(1 — 8/2) J, (1 + 8/2) J] and (J A ) = 
J(T) = $qs/167t 3 A 2 (T) with A(T) the mean-field penetration depth. We consider a field of 
/ = 1/6 and a system of size L 3 . This model was studied in Ref. ||, where it was shown that 
the VL melting transition is first order. It also has a simultaneous loss of superconducting 
coherence both in the xy planes and the z direction ||, as seen experimentally in YBa2Cus07 
I- 

The current / A (r) in each bond of the network is || 

W = *° dA f, (r) + Im sin [A A 0(r) - A A (r)] + ^(r, t) (2) 
ztiKn at 

with Jo,a = 2vrJ^/$o, and -R^v the shunt resistance. The thermal noise term is taken to have 
correlations (^(r, t)^/(r / , t')) = {2ksT J Rj l ;)5p, } [i,i8 T ^8{t — t'). Together with the condition 
of current conservation, 

£[I A (r) - I A (r - p,)] = A A ■ J A (r) = I ext (v), (3) 
A 

this determines the full set of dynamical equations. The boundary conditions are periodic 
along the x (si) and z directions, and open in the §2 direction, with a current bias / 
corresponding to I ex t( r ) = I(S r2 ,o ~ ^t 2 ,l)- We simulate this set of equations with the same 
numerical methods and integration parameters as in Ref. ||. 

We calculate the normalized voltage drop along the direction of the current as v — 
jj Sn, ,7-3 (^( r i> A r 3) ~~ $( r i, 0, r 3 )^. In Fig. 1(a) we show the dc resistance R = v/i for a high 
current i = 0.2 (currents are normalized by (io,/i)) as a function of increasing and decreasing 
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temperature, for systems with disorder 5 = 0.01,0.1,0.8. The driving current is well above 
the critical current i c (T = 0) ~ 0.04, so the moving lattice is depinned. The local vorticity 
along the direction v is given by rii>(R, t) — — Spiaquette mn t [(A/i6>(r) ~~ ^W 1 "))^ 71 "] for v _L fi; 
with nintfa;] the nearest integer to x. Vorticity is conserved both locally: A M • n^(R) = 0, 
and globally: (n z (R,t)) = f, (ng t (R, t)) = 0. In order to study vortex lattice melting, we 
calculate the vortex structure function: S(k) = J2r 3 ( n z(k, r^)n z {— k, r 3 )). In Fig. 1(b) we 
show the intensity of one of the Bragg peaks of the VL, Sq = S'(Gi), with Gi = (^, as 
a function of T for 5 = 0.01, 0.1, 0.8. Besides the field-induced vortex lines, it is also possible 
to have thermally induced vortex loops. Therefore, we calculate the average number of extra 
vortices along the z direction, n z = (\n z (R, t)\) — f, and the average vortex excitations in 
the xy directions n xy = (\n xy (R, t)\). This is shown in Fig. 1(c) as a function of T for 5 = 0.1 
only. 

The results of Fig. 1(b) show that the moving VL has a first-order melting transition 
where Sq vanishes with a sudden drop of nearly two orders of magnitude at a melting 
temperature T[j. We also see hysteresis when increasing and decreasing T. We find that for 
weak disorder 5 < 5 C ~ 0.5 the melting transition does not show any significant variation 
with disorder (in Fig. 1(b) the results for S = 0.01 and 5 = 0.1 are very similar). This means 
that the effect of random pinning in the moving VL is negligible for small disorder. For 
5 > 5 C , Tlj shifts to lower temperatures and the hysteresis loop slightly increases (shown for 
5 = 0.8 in Fig.l). A finite size study of the transition is shown in Fig. 2 for 5 = 0.1, for 



system sizes L = 12, 18, 24 ||20|| . We see that for L = 12, Sq vs. T shows a broad transition 
with little hysteresis. After increasing the system size, there is a clear hysteretic jump in 
Sq, which becomes sharper when going from L = 18 to L = 24, suggesting that there is 
a first-order transition in the thermodynamic limit. The melting temperature T[,j slightly 
decreases when going from L = 18 to L = 24, meaning that the asymptotic behavior has 
not been reached yet. (However, the difference is close to the temperature resolution of this 
Langevin simulation, AT « 0.01). We also find that the superconducting coherence along 
the z direction vanishes simultaneously with the melting transition, from the calculation of 



the helicity modulus T z (not shown here), similar to what we found in equilibrium for 5 = 



The experiments of Ref. |13| , p^] consisted of measurements of voltage and differential 
resistance for the VL driven at / ^> I c . A small dip in the differential resistance was 
observed at temperatures close to but below the equilibrium T M . In Fig. 1(a) we see that 
the dc voltage for T < T T M is approximately constant; it has a small dip at T T M and then 
rises with temperature for T > T\ t . The dip in the voltage is very small for weak disorder 
(within the statistical error), and increases with disorder, being very noticeable for 5 > 5 C . 
The voltage dip is usually associated with a peak effect in the critical current happening 
right below melting |H||14]]. From our simulation results, we see that this is only a weak 
evidence for a melting transition. Recent magnetization measurements done by Fendrich et 



al Jl5| give a better signature of the melting of the moving VL. 

It is interesting to study the behavior of thermally excited vortices. In Fig. 1(c) we show 
the number of vortex excitations in the directions parallel to the field, n z , parallel to the 
current, n y , and perpendicular to both the field and the current, n x . We see that all of them 
increase with T (at low T it is a thermally activated process n u ~ exp(—U u /T)), and they 
have a sudden jump and hysteresis at T/ f . We also find that n x > n y ~^> n z for all T. When 
compared with their values in equilibrium we find that n x (I) > n x (0), n y (I) < n y (0) and 
n z (I) > n z (0) (near equilibrium n x (0) = n y (0) ^> n 2 (0)). This means that the effect of the 
drive is to orient the vortex loops in the plane perpendicular to the current, and to increase 
the average size of the loops. Also the relative fraction of vortex excitations in the field 
direction n z /(n x + n y ) increases with J, thus increasing the entropy of the moving phase. 
These effects tend to lower the temperature for loop unbinding with respect to the melting 
and unbinding temperature of the static VL. 

Let us now analyze the structure of the moving VL. Surface plots of S(k) are shown in 
Figs.3(a-b) for temperatures below and above T l M . We find that below T[j there are Bragg 
peaks which suddenly disappear at T^. The Bragg peak structure of the low temperature 
moving phase is very anisotropic as can be seen in Fig. 3 (a). There are two possible non- 



degenerate orientations of the VL relative to the external current. In the 2D simulations of 



Moon et al. |L8] the VL is oriented with one of the reciprocal lattice vectors perpendicular 
to the driving Lorentz force, while here one is parallel to it. This orientation is favored here 
because of the network discretization relative to the direction of the applied current. (In 
the other orientation the VL is frustrated.) In Fig. 3(a) we see that the Bragg peak with 
reciprocal lattice vector Go = ±(-^,0), parallel to the direction of motion, is considerably 
smaller than the peaks at Gi = ±(^f-, and G2 = ±(— ^=). A similar result was found 
in [18[] for the peaks with components parallel to the direction of motion. However, when 



increasing system sizes, we do not see here any significant decrease on the height of the Go 



peak, as reported in [18]. In Fig.3(c-d) we analyze in detail the finite size scaling behavior 
of the peak at G2. The peak is very anisotropic: it has a finite width in the direction 
perpendicular to vortex motion [Fig. 3(c)] whereas it has almost zero width in the direction 
parallel to vortex motion [Fig.3(d)]. In Fig.3(c) we show a finite size scaling analysis of 
the peak width along the y direction : S(5k y , L) ~ L 2 ~ u F(5k y L) for L = 12, 18, 24, with 
Sky = k y — Giy- Due to the discreteness of the network there are few points to consider, 
and the scaling results are not very precise. We obtain v ~ 0.15 ± 0.1. We find a similar 
result for the other peaks, both for Gi and for Go, but with a greater inaccuracy in the 
latter case. Despite the large uncertainty, from our results we can be sure that v < 0.3. In 



the 2D simulations of [18[ an exponent of u 2 d ~ 0.53 was obtained, with seemingly isotropic 
Bragg peaks. The anisotropic peaks found here are more consistent with the discussion of 
Giamarchi and Le Doussal in [|16|], where it was shown that the effects of the static disorder 
are still relevant in the direction transverse to vortex motion. The fact that in our case 
the low temperature moving phase seems more "ordered" (i.e. smaller v) can be due to: 
(a) three-dimensionality; (b) the long-range (logarithmic) interactions between vortices. In 
any case, our results are consistent with the recent experiment of Fendrich et al. []nj where 
a large external drive does not affect the first-order character of the melting transition, 
indicating that there is a low temperature moving solid. 

We have made a systematic comparison of our results for the moving VL (I 3> I c ) 
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with simulations near equilibrium (J <C J c ) for different values of 5. The main result is 
summarized in the phase diagram shown in Fig. 4. We see that for increasing disorder the 
equilibrium first-order transition turns into a continuous transition or a crossover for 5 ^ 0.5, 
in agreement with the results of Jagla and Balseiro ||. On the contrary, the phase transition 
of the moving VL is always of first-order character, the only effect of large disorder 5 > 0.5 is 
in lowering T M . Also we find that T\ { < T M for any value of the disorder. This is due to the 
coincidence of melting with vortex loop unbinding: since the bias current tends to unbind 
loops more easily, T M should be lower than T M , as we discussed above. Also at zero field, 
where the phase transition is driven by vortex loop unbinding, a bias current considerably 
lowers the critical temperature plf . Recently, it has been suggested by Nguyen et al. 



that the presence of vortex loops in the melting transition could explain the anomalous 
behavior of the entropy jump found by Zeldov et al. [|]]. From our results, Fig.4, we believe 
that a measurement of a finite ATm = T M — T M could be an experimental probe of the role 
of vortex loops in the melting transition. Also a "shaking temperature" effect could give 
T M <T M JT2|. However, this gives a ATm which depends on disorder, and should vanish 
for 8 — > 0. In the experiment of Safar et al. [13], with a field B || c, a dip in the differential 



resistance was observed very close to, but below, the equilibrium melting transition. The 
experiment of Fendrich et al. []U| does not show a reasonable dependence of Tm with the bias 
current. This can be due to: (a) the ATm is too small to be discernible within experimental 



resolution, or (b) the bias current of [|TfJ is not high enough (here, the current is such that 



the VL is in a flux flow state at T = 0, whereas in |l5j there is no dissipation at low T). 



In the experiment of D'Anna et al. ||14j| , with B || ah a dip in the differential resistance 



appears clearly below the equilibrium Tm- Since it is easier to have thermal nucleation of 
loops parallel to the ah planes, this could lead to a measurable ATm for B || ah. 

In conclusion, we find that the low temperature moving phase is an anisotropic glass 
which melts via a first-order transition in 3D, even for large values of disorder. Our re- 
sults stress the three-dimensionality of the VL: there is a more "ordered" moving phase 
when compared with 2D simulations fl8[], and at melting there is a simultaneous loss of 



superconducting coherence along the z axis. 

D. D. acknowledges discussions with H. Safar. This work was supported by the 
U.S.D.O.E.. 
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FIGURES 

FIG. 1. (a) DC voltage over current R = v/i as a function of temperature T for a bias current 
i = 0.2 3> i c and L = 18. Dotted line: 5 = 0.01, triangles: 8 = 0.1; squares: 6 = 0.8. The melting 
temperatures in equilibrium, T^, and for the driven system, T[f, are indicated, (b) Intensity of the 
Bragg peak Sq vs. T. (c) Plot of the density of thermally excited vortices n v vs. T, for S = 0.1. 

FIG. 2. Plot of the intensity of the Bragg peak Sq vs. T for S = 0.1 and different system sizes. 

FIG. 3. Surface plots of the structure factor S(k) for S = 0.1, current i = 0.2 and L = 18. For 
(a) T = 1.0, (b) T = 1.2. (c) Finite size scaling plot of the k-dependence of the Bragg peak in 
the direction transverse to vortex motion, (d) Same as (c) but in the direction parallel to vortex 
motion. Results for S = 0.1 and T = 0.9. 

FIG. 4. Disorder-temperature phase diagram. Triangles: equilibrium transition temperature 
Squares: transition temperature of the driven system T^. Continuous lines: first-order 
transitions. Dashed line: continuous transition. Dotted lines: limit of the hysteresis loops for a 
given cooling rate. Results for L = 18. 
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